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Abstract 

We show how kinetic theory, the statistics of classical particles obeying Newtonian dynamics, 
can be formulated as a field theory. The field theory can be organized to produce a self-consistent 
perturbation theory expansion in an effective interaction potential. The need for a self-consistent 
approach is suggested by our interest in investigating ergodic-nonergodic transitions in dense fluids. 
The formal structure we develop has been implemented in detail for the simpler case of Smolu- 
chowski dynamics. One aspect of the approach is the identification of a core problem spanned by 
the variables p the number density and B a response density. In this paper we set up the pertur- 
bation theory expansion with explicit development at zeroth and first order. We also determine all 
of the cumulants in the noninteracting limit among the core variables p and B. 







I. INTRODUCTION 



There exists a well defined approach to the problem of classical many-particle dynamics. 
Kinetic theory governs the kinetics of particles obeying Newtonian dynamics. It is one of the 
oldest disciplines in all of science. From the early work of Bernoulli^ to the seminal work 
of Boltzmann2] and Maxwelljsl, kinetic theory has been applied to dilute systems out of 
equilibrium. More systematic modern methods have been developed for systems fluctuating 



in equilibrium. Based on work by Koopman 



J] and von Neumann 



5|, one can express the 



time evolution in terms of Koopman's operator e*^* where L is the Liouville operator. For 
example one can develop in a rather straight forward way density expansions for transport 



coefficients 



and memory 



of the interaction potential 



unctions 



Similarly one can develop expansions in terms 



10|. However these approaches have their limitations. It is 



not known how to systematically rearrange the respective density or potential expansions 
in a self-consistent manner. By self-consistent we mean here the interaction kernels of the 
associated kinetic equations can be expressed in terms of the full unknown correlation func- 
tions. Thus one obtains a nonlinear kinetic equation that must be solved self-consistently. 
While there are many formal reasons why this self-consistency is desirable, our motivation 
for pursuing such a theory in the present case is more practical. We want to understand 
the role of ergodic-nonergodic (ENE) transitions fll|. Self-consistency is essential if one is 
to investigate whether one does or does not have an ENE transition in dense fluids. 

The liquid cooled to low temperature or compressed to high densities reach a stage 
in which it behaves like a frozen solid without any long range order. Understanding the 
formation of the amorphous solid state of the liquid coming from the ergodic liquid state 
has remained an unresolved problem of physics. A theory of the formation of the amorphous 
solid like state with self generated disorder will require developing techniques for testing the 
possibilities of ENE transitions in the dense liquid. There is compelling evidence that as one 
approaches the glass transition one comes close to an ENE transition. Even if we finally rule 
out the existence of a physical ENE transition, we need a formal structure which potentially 
can show such a transition. We present the outline of such a self-consistent theory here for 
the case of Newtonian dynamics (ND). Elsewhere one of us has developed 



12 



13j such a 



self-consistent theory for a classical set of particles which obey the simpler Smoulochowski 
dynamics fl3|(SD). In that case a program of investigating the existence of ENE transitions 
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is well along and one has a well defined approach to the problem. A key attribute of this 
theory is that well-defined approximations are available with well-defined corrections. In 
this paper we show at the formal level that Newtonian dynamics can be organized in the 
same self-consistent fashion. Our main goal is to investigate the status of ergodic-nonergodic 
transitions in dense liquids. How does one approach the liquid-glass transition? In Ref. 
the question of an ENE transition in a system obeying SD by working at one-loop order was 
addressed. With the simplest interaction vertices it was demonstrated that a system of hard 
spheres does not under go an ENE transition until a packing fraction of rj* = 0.76. This is 
the reult of the simplest calculation. There are a variety of perturbations and extensions 
one can workout to test the robustness of this result. A better description for the static 
structure and/or improvements of the vertex functions to include interaction effects can also 
be used in the calculation. Going to higher loop order is also another way of improving 
the predictions of the theory. We anticipate that all of these calculations will eventually be 
carried out for SD. 

There are differences between Newtonian and Smolukowsi dynamics. First, Newtonian 
dynamics are reversible and Smoluchowski dynamics have a noise component. Related is 
the fact that ND has additional conservation laws (energy and momentum) compared to 
SD. Finally there is the technically very important fact that the phase-space in the ND case 
contains the particle momenta while one has only the positions in the SD case. In the SD 
case the kinetic equations lack the momentum index labeling labeling ND kinetic equations. 
This is a large technical advantage of working with the SD. In the present work we take 
the point of view, which can be checked, that as one approaches the glass transition the 
slow kinetics are associated with structural rearrangements of the density, and couplings 
to energy and momentum currents are less important. This is one of the assumptions of 
simplest version of the mode-coupling theory [15j]. A related working assumption is that 
the glassy kinetics of SD and ND are very similar. We know that they share the same 
static correlation functions-static structure factors. We will eventually be able to show this 
similarity. 

An important formal point is that within the theory there exists a core problem involving 
p, the particle density, and B, a response density. One must address this core problem in 
both SD and ND before including other variables like momentum variables in ND. More 
specifically a discussion of the rich hydrodynamical structure of this system is addressed as 
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one introduces additional variables. For example the shear viscosity is associated with the 
inclusion of transverse currents in the development. 

In the next section we show how kinetic theory can be reexpressed in terms of path 
integrals. The path integral formulation is similar to that found in the SD case with the 
initial condition playing a role similar to that for the noise in the SD case. In the SD case 12| 
the path integral approach involved the same core field s: p and B. The ND case can also be 



organized most simply in terms the same two variables [16|]. With respect to their respective 
interaction structures, the formulations are same for ND and SD dynamics. We find that 
the noninteracting cumulants for the set of variables p and B are qualitatively similar for the 
two cases. Thus we follow here the development in Ref. 12| and construct a self- consistent 
approach to the dynamics. 

It has turned out in the SD case that we can replace the bare potential with an effective 
potential expressed in terms of the physical structure factor. Thus the theory can be applied 
to systems with hard-sphere interactions. Here we develop a perturbation expansion for the 
two-point cumulant in terms of the pseudo interaction potential. The glue which ties together 
the terms in perturbation theory are the three-point vertices which are constructed from the 
noninteracting three-point cumulants among the variables p and B. We show that all of 
these cumulants can be evaluated in the time and wave number regime. We focus here on 
the two-point cumulants and their determination to first order in the interaction including 
the renormalization of the bare potential. 

At the formal level we show that the collective contribution to the self-energies for p 
and B have the same form in the ND case as for the SD case to second order in the 
expansion. The main calculation carried out in this paper is the determination of all the 
cumulants between p and B for the noninteracting system. In the next paper in this sequence 
we address the development of the theory at second-order in the interaction. The first 
important development there is to show, more simply than in the SD case, that one has a 
fluctuation dissipation symmetry relating the response functions to the correlation functions. 
In this case one can establish a set of nonperturbative identities satisfied by the three-point 
cumulants and irreducible vertices. Additional identities are clearly available for higher 
order quantities. In turn one shows that the set of matrix Dyson's equations, satisfied by 
the two-point cumulants, reduce to a single kinetic equation for the correlation function as 
assumed in mode- coupling theory. It is from this kinetic equation that one can develop the 
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machinery associated with the ENE transition and the slow dynamics one can observe even 
if one can not access a sharp ENE transition. 

In this paper we lay the ground work for the analogous calculations for Newtonian dy- 
namics. These calculations are in some ways harder in the ND case. For example, the 
three-point vertex functions are more complicated. However, in other ways the evaluations 

he expansion param- 



13| mentioned above. 



are easier. The treatment of the FD symmetry is simpler in ND case, 
eter in this problem is a pseudo-interaction potential just as in Ref. 
Indeed we expect to find the same interaction pseudo-potential in the two cases. 

At second order we have an expression for the static structure factor S{q) in terms of 
the potential. We assume we can pick the best result for S{q) and solve for the associated 
pseudo-potential. Plots of V{q) for various packing fractions ( e.g., see eqn. (218) and Fig. 2 
in Ref. [l3| ) shows that the pseudo-potential is a rather smooth quantity even for rigid 
hard spheres potentials. 

A key difference between SD and ND is the noise in the SD case driving the system 
toward equilibrium. At low orders in the ND system it needs to be told it is in equilibrium. 
In conventional kinetic theory the initial state is chosen to be in equilibrium. Once in 
equilibrium the system typically stays there. Thus one option is to fix an initial condition at 
t = to. We show how this works in the zeroth and first-order cases for the dynamic structure 
factor. In the first-order case one has two effects on the theoretical structure. First one has 
(see Eq.( !27l) below) that the interaction matrix has a piece which imposes an initial condition 
and, second, the time integrations are restricted to the time regime to < t. We show that 
the fluctuation-dissipation theorem[l^] (FDT) is crucial in seeing that these breaks in time- 
translational invariance (TTI) cancel out and one obtains correlation functions compatible 
with TTI. In our treatment of FDT we flnd that the system is invariant over a set of 
symmetry operations which depends on an undetermined parameter. We propose to flx 
the system at temperature T by choosing that all cumulants and vertices satisfy the FDS 
associated with equilibrium. 

II. NEWTONIAN DYNAMICS 

Consider a system of particles with mass m with conflgurations specifled by the phase- 
space coordinates = {Ri,Pi) which satisfy the equations of motion: 
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where the particles experience force 



f^ 



_d_ 



U{R) 



(3) 



with total potential energy 



U[R) = -Y.V{R,-Ri, 

i+3 



(4) 



and we have suppressed vector labels to unclutter the equations. If we form the vectors 
= {Ri, Pi) then the equations of motion ([1]) and ([2]) can be put into the form 



(5) 



where i^j is a : 
MSR structure 



unction of the \E'. We treat the phase-space coordinates as our fields in a 



18| and the physical observables are treated as conjugate to external fields 



which label the generators of the physical observables. The A^-particle partition function is 
given by 



Z^[H,h,h] = M \{V{mi)V{i>,)d^f'^)Po{^,)e 

i=l 

X exp(iJ ■ 0) exp(/i ■ \E' + i/i • ^) 



(6) 



where A/" is a normalization constant, and we have an initial probability distribution Po(^o)- 
We assume the system is in equilibrium initially and the initial distribution is canonical: 



(7) 



where "H is the hamiltonian, which is the sum of the kinetic and potential energies. The 
MSR action for the problem is given by 



Aq, = dt 



(8) 



Finally we have the contribution in Eq.® due to external fields that couple to $ the 
collective core variables of interest: 



H-^ = Y^j dlH^{l)^^{l) . (9) 



As explained in Ref. [12|, the minimal set for $ includes the particle density and the 
response field B. For convenience we have also included the fundamental source fields hi. It 
is useful to show that this development can be mapped onto the traditional representation. 
Set if = in Eq.([H]). No information is lost as long as we keep the full set of external fields 
h: 

Z^h] =U j P(^)P(^)d*(°))Po(*o)e['*(^+'^)] exp(/i ■ ^) (10) 

where in each of the arguments of the exponentials on the RHS there is a summation over 
particle label, an index labeling position and momentum, and an integral over time. Now 
we do the functional integral over ^ to obtain a functional 5-function: 

ZN[h] = M j I?(^)c/^(°))Po(*o)^ + k) exp(/i ■ ^) . (11) 

The next step is to recognize that in a deterministic system the probability of finding the 
system in configuration \l/(t) after starting at time to in configuration \l/o is proportional to 

5(^(t)-^(t;to)) (12) 

where \l/(t;to) is the unique configuration at time t evolving from ^(to)- In the operator 
formulation 

^(t;to)=e*''(*-*°)^(to) (13) 

where C is the Liouville operator. As discussed in some detail by Penco and Maurojiol. one 
can use the following argument to connect our development here to the standard formulation. 
If a function /(0) has a zero at = 0o then 

5(0-0o)=5[/(0)]|/'(0o)| (14) 
In our case we choose / to be \I/ + /T, we have the identity 
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S - ^{t; to)) = 6 (^(t) + K{t)) AT (15) 
and the factor A/" in Eqs.f llll) and f|T5l) is the functional determinant 



5[^{t) + K{t), 

This quantity is essentially the Jacobian discussed in Ref. [12)] which opens up the discussion 
again to ghost fermions, supersymmetry, unification [2^. Importantly, for our purposes here, 
A/" is a constant 2l| independent of \E'. Using Eq. (fT5|) in Eq.f llip leads to the partition 



function as 



ZN[h] = /l^(v|/)(iv|/(o)Po(vI/o)[]5(vl/(t)-v]>(t;to))exp(/i-v|/) . (17) 
We can immediately do the functional integral over \E'(t) to obtain 



N 



ZN[h]= j d¥''^Po{^o)exp{r dtY,hi{t)-^>i{t-to)) . (18) 

Clearly by taking functional derivatives we can generate the average of any set of phase-space 
observables: 



{A{t^)Bit2) . . . D{tn)) = I d¥'^Poi-^o)Ai-^{h; to))5(^(t2; to)) • • • ^(^(^n; to)) (19) 

and we see that our representation is equivalent to the standard theory. 

The appropriate generating functional for the problems discussed here, working in the grand 

canonical ensemble, is given by 

W[H, h, h] = InZriH, h, h] (20) 

where 

ZT[H,h,h]=Y.(^Z4H,h,h] (21) 

and Zn[H, h,h] is given by Eq.flTSll. 

For the method to be effective we have a minimum of two collective fields dictated by the 
structure of the interactions. One essential field is the particle density 
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N 



pil) = J25{x.-R.ih)) (22) 



i=l 

and it is crucial to include the response field 

TV 

5(l)=^[(A(ti))^Vi]<5(xi-i?,(ti)) . (23) 
1=1 

Notice that B depends on the MSR hatted field Pi{t). Unlike the SD case the Jacobian does 
not contribute to the definition of B{1) here. 

We can then write the canonical partition function given by Eq.dH]) in the form 

Zn[H, h,h]= H V{^,)V{^i)V{^f^)e-^''-^'+^-'^+^-^+^-^ (24) 

i=l 

where Aq is the quadratic part of the action including the quadratic contribution to the 
initial probability distribution 



N 

(0) 



^0= / dti ^i^4- ['^^ + K 



+ l3Ko . (25) 



Kq is the initial kinetic energy. Notice that we have constructed things such that the phase- 
space variables are constrained to their initial values at t = to- We then average over these 
values. Here we are explicitly treating the case where the system is in equilibrium at t = to, 
but more general situations are clearly compatible with the development. The interaction 
part of the action is given in the compact form 

^i = IT. [ dld2^^{l)ac,,{12)^,{2) (26) 

where the Greek labels range over p and B and the interaction matrix is defined just as in 
SD,by 

a^p{l2) = V{12) %p5ppP5{ti - to) + d^sS^p + K^^b] (27) 

where the first contribution is from the potential energy contribution to the initial condition 
and 

V{12) = V{xi-X2)6{h-t2) . (28) 
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Notice that the response field B is chosen such that the interaction part of the action has 
the form given by Eq.f l26p . The canonical partition function can be written in the 
convenient form 



Zj^ = fre-^'+^-'^ (29) 
where we have introduced the average 

r ^ 

Tr = / []r'(^i)r'(\^i)ci^f e-^« • (30) 

i=l 

We have thus shown that ND can be written in a path-integral form as given in 
Eqns. (]B15|) - (]B17p . These expressions look like the formulation for SD with an important 
exception. The action A^, given by Eq.®, does not have a contribution from noise. A noise 
component, as in Fokker-Planck dynamics, adds a term to the action 

roo roo ^ ^ 

Anoise= dh dt2Y.^i{h)D{hM)Ait2) (31) 

J to J to 

which contributes to the noninteracting part of the action Eq.(2.20). The noise has the 
property of continuously telling the system to equilibrate at temperature T. How in the 
case of Newtonian reversible dynamics does the system know it is in equilibrium? One 
mechanism is to satisfy an initial condition. In the conventional formulation an equilibrium 
correlation function is given by 

CAB{t) = TrPoBe^'^'A (32) 

where C is the Liouville operator and Pq the equilibrium probability distribution 

and at t = one has explicitly an equilibrium probability distribution. The use of initial 

conditions is one important way of treating nonequilibrium kinetics as discussed in Ref. 



12| . The use of initial conditions in the case where the system is in equilibrium for all 
times, including the t = Iq, is inconvenient since formally it looks like one has broken time- 
translational invariance (TTI). As shown below for an ideal gas and to first order in the 
interaction, one can tell the gas its in equilibrium with an initial condition in equilibrium 
and maintaining TTI. However there is a cleaner way of maintaining equilibrium and TTI. 
This is to require that the fiuctuation-dissipation symmetry hold at all times. 
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A. Fluctuation-Dissipation relations 

The correlation functions, via its definition as the average of the product of commuting 
(classical) fields satisfies the time reversal symmetry 

Gij{t - t') = Gji{t' - t) (33) 

The physical fields are real so G*j{t — t') — {'^i{t)'^j{t'))* — Gij{t — t'). In a stationary state 
we have time translational invariance. The full MSR action as out hned above is obtained 
in the following form for the Newtonian Dynamics 



Mt){Riit) - ^} + im{m - Fi{t)} 

TJX 



dt 



(34) 



Under complex conjugation the MSR action transforms as 



We introduce the following transformation: 



(35) 



TR,{t) = Ri{-t) 

TRi{t) = -Ri{-t)+i/3Fi{-t) 

rPiit) = -A(-t) - -P^{-t) (36) 

where Fi{t) — {dU / dRi{t)) is the force on the particle i. We now consider how the action 
A changes under this transformation 



X = tA 



iRi{-t) 



dRii-t) Pjj-t) 
d(-t) m 



+ 



+°° dH{-t) 



d{t) 
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Looking at individual terms, and letting f — > — i in the integrals, we obtain A' — tA — 
A — /3H{— oo)+ /3H {+00). We take the two limits of time integration ^2 — >■ 00 and ti — >■ —00 
here. Treating the last part as a constant (=0 in a conserved case) we conclude that the 
MSR action remains invariant under this transformation. 



B. FDT involving the p and B fields 

Let us consider the transformation rule for the field B{x,t) under time reversal r 

B{x, t)^-Y.P,- -^6{x - R,{t)) (37) 
i=i oHi 

In the present formulation we are working so far with the density variable t) as the only 
collective variable and B{x,t) is the hatted counterpart. Under r the field B{x,t) changes 
as 



TB{x,t) = -B{x,-t)--(5^^^^'''^^ 

m 



d{t) 

Using this relation we obtain for any function f\p\ the following FDT relation hnking to the 

B{x,t) 



GMt-t') = ^e{t-t')p^Gf,{t-t') 



(38) 



III. SELF-CONSISTENT DEVELOPMENT 



The self consistent theory for the liquid state dynamics is developed in terms of correlation 
functions of collective variables S introduced in the previous section. We use the notation 
for the variable as 

i 

with the index a denoting the space of collective variables and i is the particle label. We 
work here with the density variable p and conjugate variable B{x, t) which are respectively 
denoted as: 
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and 



(40) 



(/.f(l)=zA(ti)V.,(5(xi-i?,(ti)) . (41) 

Working in the grand canonical ensemble, the grand partition function for the interacting 
problem is given by Eq. fl2T]) . The cumulants of the fields $j are generated by taking 
functional derivatives of the generating functional 

W\E\ =ln Zt . (42) 

with respect the fields H introduced in eqn. ([9]) above. The one-point average in a field is 
given by 

G. = (43) 

In the above equation we have used a compact notation where the index i labels space, time 



and fields p or B. We maintain the notation from here on. In Ref. 12|] we derived the 
fundamental identity 

Gi = Tr0,e^-<^+^^ (44) 

where 

AW = W[H + F]- W[H] (45) 
with H denoting the external field and Fi being given by 

Fi = Y.ai^ . (46) 

j 

with the interaction kernel o"jj obtained in Eq. f l27|) . The important result Eq. fj44|) was 
established in FTSPD. It is more useful to derive here this result in a completely different 
fashion. This is presented in Appendix|A]in which we derive Eq.f H^ without using functional 
techniques. 
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IV. FORMAL DEVELOPMENT OF PERTURBATION THEORY 

The dependence of the theory on the pair potential is controlled by the quantity 



Aiy = W[H + F]- W[H] (47) 

in eqn. (I44p where F is proportional to the interaction potential. We can expose the depen- 
dence on the potential by constructing the functional Taylor-series expansion 

AW' = E F.^Wim + J: If.F.^WIH] + . . . (48) 

and we can introduce the full cumulants: 

Gii k = — — —...^W[H] (49) 

to obtain 

AW = Y^ F^G, + ^F,F,G,, + ^ ^F.FjFkCj, + ... . (50) 

i ij ijk 

Clearly in this form we can take AW to be a functional of Gi. One can then use functional 
differentiation to express higher order cumulants in terms of Gi and Gij . One has for 
example the manipulation expressing the three-point cumulant in terms of Gij 



Gijk r- jj Gij ^ ' Gifj^G jnG]^p^ Ynnp (^'^) 
"^k mnp 



and the irreducible three-point vertex. Tij^ in turn is given as a functional derivative of the 
two-point irreducible vertex, 

^ijk = -p^^ij (52) 
^k 

with respect to Gk- ^ij is precisely the matrix inverse of the two-point cumulant Gkj'- 



^ikGkj — Sij (53) 

k 



where we refer to this as Dyson's equation. 



13 



As in the SD case 13|] we can establish a dynamic generahzation of the static Ornstein- 
Zernike relation[22|. Starting with the functional equation for the two-point cumulant, we 
use the chain-rule for functional differentiation to obtain: 



— c 



'5ij + CikGkj (54) 

k 



where 



is roughly speaking a one-body object, and Cj, is defined as 



g,j = Tr<^,</),e^-<^+^^ (55) 
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- Tr0,e^-^+^^-^AW^ . (56) 

Since AW can be treated as a functional of Gi we see at this stage that we have available 
a self-consistent theory. To formulate this we define the matrix-inverse of the one body 
quantity Qij as 'jij, 

X! likGkj = Sij . (57) 

k 

Now multiplying the correlation matrix G defined in eqn. flS^ with the matrix 7, and using 
definitions fl53l) and fl57|) . we obtain the two-point vertex function (without any approxima- 
tion) : 

^ij = 7ij + ^ij (58) 
where the collective part or the dynamic self-energy is given by 



Kij = - ^ likCkj ■ (59) 

k 
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A. Collective Self-Energy at First and Second Order 

At the formal level we can work out the collective part of the self-energy in perturbation 
theory. Using Eq.f l56p we have at first order 

c« = rr0,e^<^+^^-|^AiyW . (60) 

where 

Aiy«=5]F,G, . (61) 

i 

Clearly 

^AW^W = F, (62) 

and 



k 

= Gik(^kj ■ (63) 



The collective contribution to the self-energy at first-order is 



Next, we consider the second-order contribution 



Kij^ = -lieGek(^kj = -(^ij ■ (64) 



where 



It is shown in Ref. 



eg) = Tr0,e^-^+^^-|^AH^(2) (gg) 



Aiy(2) = . (66) 

13| that the second order contribution to the dynamic self energy is 



obtained in the symmetric form in terms of the screened matrix correlation function G as 
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where '-yf-]. denotes the three point vertex function Tijk defined in eqn. ( I52l) at the lowest order 
i.e., for a noninteracting system. This is referred to as the collective one-loop contribution. 
The G matrix elements are obtained as 



GijiX) ~ 9 X! \pfk(^kiGij + GikCTkiG^f (68) 



^ kl 



and plays the role of an effective propagator. The above result obtained for the SD case 



in Ref. [l3j also holds for the ND case. The relations ( 158|) and ( 167|) form a closed set 
of self-consistent equations for the correlation matrix Gij. This gives rise to a feed back 
mechanism which becomes strong enough at high density to drive the system to a possible 
ENE transition. Evaluating the strength of the feed back mechanism even at the lowest 
order will involve evaluating the three- point vertex functions ^y^jl for the noninteracting 
system. Later we describe the basic calculation for obtaining the correlations in the ideal 
gas. Analysis of the feed back process with evaluation of the vertex functions and the 
possibility of the ENE transition will be taken up in a companion paper. 



B. First-Order Theory 

To obtain the correlation functions at the first order we need as an input the correlation 
functions for the noninteracting system. We present computation of the two point correlation 
functions G^^j at the zeroth order in the next section. We list the zeroth- order results here 
starting point of the first- order calculation : 



^^(^,^1,^2) 


= Poe 2/3'" 






(69) 




k\ 

= ti 

m 




- t2)p0 


(70) 




k\ 

= — ih 

m 




-ti)po 


(71) 




= 






(72) 



where we have the same structure as in the SD case. Gpp is real, GpB is retarded, Gsp is 
advanced, and is zero. For the time translational invariance we obtain the frequency 



transformed quantities as 



,2 



GfJ{k,u) = V2^^e ^ (73) 
G^Sik,.) ^ (74) 

where ttiVq = ksT and where the integral S{x) is defined as 

S{x) = 1 - 2x'e-"' £ due""' + i^xe"'^ (75) 

Let us look at the first-order theory for the one point function or the equation of state. 
Starting from the basic relation (144|) we obtain for the one point function 



= Gf)+G;°VHG, (76) 

For i = p using the corresponding a^p element and using the result Gp{l) = 2TT6{ui)Gp{ki) 
we obtain 



Gpik) = ^ (77) 



l-Vik)G%{k,0) 

Since G^p\0) = po we obtain the equation of state to first order in V as 



where we have used the result from eqn. (fMll that G||^(0, 0) = — /3po to zeroth order in the 
interaction. Notice that we self-consistently determine the static properties as we solve the 
dynamic problem. 

It turns out to be more useful to write eqn. (1781) in the form: 

^° = 1 - PpV{0) (^'^ 
= pexp{V{0)) . (80) 
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We now relate this to the conventional equation of state. First, since we are working in the 
grand canonical ensemble, po is the fugacity 

Po = ^ (81) 

where £ is a microscopic length independent of density and p is the chemical potential. There 
is the Gibbs-Duhem relation 

or 



dl3p _d(3fi 
P dpo 



Po dp 

Working at first order in the potential it is easy to show 



(84) 



dpo V 



dp 

Then 



e^loj(l + y[0]) (85) 



^ = (1 + V[0]) (86) 

/3p = P+^/3pV[0] (87) 

The first term is the ideal gas law. The complete expression with V an effective potential is 
sensible over a broad range of densities. 

Let us look at the first-order theory for the two-point correlation functions. To obtain the 
correlation functions at this order we begin with the eqn. ( |5^ and use the corresponding 
order expression for the self energy cfj^ and AM/*^^) obtained respectively in eqns. ( 163|) and 



W\\\ to write the matrix kinetic equation 12|] : 



Gij - G\f + cf^iCTikGk + X! G'^a^^ikGkj (88) 



The role of the second term on the RHS of the above equation gives rise to a factor 
p/ Po with the first term. In time space we write. 



12| of 
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Gi,{k,h,t2) = ^Glf{k,h,t2)+ r diY,Gfl\k,h,t)aM{k)Gi,{k,i,t2) . (89) 

The average density p is defined in the equation of state given by eqn. ( 1781) in the time 
translational invariant form. Now from eqn. fl89p we obtain setting i = j = B and t2 = to 
(the initial time) 



GBBik,h,t2) = -G%{k,hM) + G%l{k,t,,to){-(3V{k))GpB{k,toM) 
Po 

(90) 

Since G'^^(A;, ti, ^2), G^l^p{k,ti,to) and G^p^{k,to,t2) are zero, we find GBB{k,ti,t2) = 0. 
Consider the response functions GpB{k,ti,t2) and GBp{k,ti,t2) by taking the p — B and 
B — p matrix elements of Eq. (l89|) results in the equations 

G,B{k,h,t2) = ^G%{k,U,t2) + r dtG%{k,U,t}V{k)G,B{k,t,t2) 

Po Jto 

GBp{kMM) = ^Gfp{kMM)+ H diG%l{kM.t)V{k)GBp{k,tM) 

Po •'to 

(91) 

It is easy to see that GpB^k, to, ^2) and GBp{k, ti, to) are both equal to zero. Since GpB^k, i, ^2) 
is retarded and G^^^p{k,ti,i) is advanced, eqn. (l9Ti) are solved by Fourier transformation: 

"''^"'^^^ l^Vik)G%ik,.) 

l-V(k)G%(k,L0) 

(92) 

Using the results from eqn. fl74|) that G'^'g(0, k) = — /3po, to zeroth order in the interaction, 
we obtain the zero-frequency GpB at the first order as. 



1 + l3poV{k) 

G'^>{k,0) = - 



^^ /5P 



1 + /3poV{k) 

(93) 
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Now we consider the p — p matrix element of Eq. (189|) 



Gpp{kMM) = ^Gf^{kMM) + Gf^{kMM){-PV{k))G,,{k,t,M) 
Po 

+ / dtG'^^J{k,h,t}V{k)GBp{k,t,t2) 

+ r dtG%{k,t,,t)V{k)G,,{k,t,t2) (94) 
We rewrite this equation in the following form 



G,,{k,h,h) = ^Gf}{k,hM)+ r dtG^^J{k,h,t)V{k)GBp{k,t,t2) 

Pq J —oo 

+ r dtG%{kM.i)V{k)Gpp{k,tM) + K{kMM) 

J —oo 



(95) 



where all of the terms dependent on are included in 



K{kMM) = GfXk,ti,to){-mk))Gpp{k,to,t 



pp 

to 



2) 



diGfJik,h,i)Vik)GBpik,i,t2 



pp 

— oo 
rto 



- (' diG%{kM,t)V{k)Gpp{k,iM) ■ (96) 

J —oo 

All of the elements contributing to the breaking TTI are collected into K. By direct 
calculation one can show that K = Q and TTI holds. This analysis is rather involved. 
Instead we consider the case where the FDT holds 

GpB{t-s) = e{t-s)l3^G,,{t-s) . (97) 

Now in K eliminate the response functions in terms of the density-density correlation func- 
tion. Notice that the two terms involving an integration can be combined into a single term 
where inside the integral one has an exact derivative. Doing the integration the contribution 
from i = to cancels the term from the initial condition and one has K = 0. Eqn.( l95l) with 
K = can be solved by Fourier transformation and one has 



G 



ip/Po)G^S 
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Vik)GPp 
1 - Vik)G^Z 



(98) 



and we have used the previously determined result for the response functions. After a little 
rearrangement we have the real result 



[^-V{k)G%\[l-V{k)GZ 

where 



(100) 



This is the solution to our first order problem. 

There are other ways of looking at the information contained in the solution. If the FDT 
holds we have in Fourier space 

^ImGpB = (101) 
which also holds for the noninteracting limit 

^ImG% = Gf) (102) 

These lead to an integral form for the FDT. We start with the identity for the Laplace 
transform of GpB 



GpB{z) 



°° duImGpB 

L 



oo IT Z — UJ 

°° du „ G 



CO- 



oo 27r z — UJ 
°° duo ^ , .UJ - z + z 



p du L 



-oo ZTT Z — OJ 

j-oo Zn z — UJ 

= -f3pS{k) + zf3Gpp{k,z) (103) 

where the Laplace transform of the density-density correlation function is defined as 

G„(.,.)=r^5^ . (104) 

J-oo ZTT Z — UJ 

We have used above the result 
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-oo ZTT 

= -l3pS{k) , (105) 



where S{k) is the static structure factor. Using eqn. ( 193|1 we obtain our first order approxi- 
mation for the structure factor 



which is equivalent to taking the two point Ornstein-Zernike direct correlation function as 



c{k) = -(3V{k) . (107) 
We then equate the two expressions for Gpsik, z) obtained in eqns. (192|) and (11031) 



GpB{,z) = -f3ppS{k) + zf3Gpp{k,z) 



{p/po)G%{k,z) 



l-V{k)Gf^{k,z) 



(108) 



Solving for Gpp we find 



Gpp{k, z) 



S{k) + 



(109) 



l-V{k)Gfi{k,z) 

In the noninteracting limit, we define ip{z) as the density autocorrelation function normalized 
with respect to its equal time value 



ij{z) = GfJ{k,z)/po . (110) 
Laplace transforming the FDT relation (11021) for the noninteracting system we obtain 



G^'!,{z) = -Ppo{l-z^{z)) 



;iii) 



Using this we obtain from eqn. (11091) the result 



Gpp{k,z) 



+ 



(1 - zi^iz)) 



1 + /3poV{k) 1 + PpoV{k){l - zij{z)) 



(112) 
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After some algebra we have for the density correlation function, 



Gpp{k,z) = pS{k) 



l + PpoV{k){l-z^{z)} 



(113) 



where 



CO 1 / ^9 

/o 



Let us look at this in the small k limit. We easily find 

(kvo)' 



i;{z) 



1 + 



(114) 



(115) 



plus terms of order /c^. Putting this into Eq. flll3p and rearranging we find 



Gpp{k, z) 



-pS{Q) 



1 1 

+ 



z — ck z + ck 



(116) 



The two poles represent the two propagating sound modes respectively with speed c{k) is 
obtained as, 



c{k) = vo{l + f3poV{k)y/' 



(3mS{k) 



(117) 



V. MOMENTUM VARIABLES, CURRENTS, AND THE HIERARCHIAL 
STRUCTURE OF THEORY 



A. Additional Degrees of Freedom 

So far the theory has been set up to deal with what we will call the core problem. This 
is the determination of the observables involving the core variables $o = {Pi B). In the case 
of SD this covers essentially all of the degrees of freedom of interest. For ND this is not the 
case. We have a number of additional degrees of freedom. This is because we have a larger 
phase-space due to the momentum degrees of freedom. Choices of variables to be included 
from the simplest to the most involved are: 

1. Couple to the two transverse currents g±,i, g±,2 

2. Couple to the whole current g and/or the kinetic energy density gx- 
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3. Couple to the phase-space density 

N 

f{x,p,t) = Y.5{p-Pi{t))S{x-Ri{t) . (118) 

1=1 

Choice 1 is the simplest since there is no direct coupling between the longitudinal and trans- 
verse degrees of freedom. This is also the simplest way of determining the shear viscosity. If 
one goes to choice 2 things are more complicated since we have all the correlation functions 
among p, B and g. One needs g,K to determine the thermal conductivity. 

If one wants to investigate the Boltzmann equation and momentum distributions one 
needs to include the phase-space density f{x,p,t) as one of the variables. This will be 
discussed elsewhere. 

B. Generic Inclusion Of Additional Variables 

Suppose that we have a single-particle additive variable Qa we want to include in our set 
$ = (p, B, g). This inclusion simply involves an additional term in the basic action 

A-H-^(^^A-H-^Q-J-g (119) 

which introduces a new external coupling into the problem. The fundamental generating 
functional is now given by 

W[H, J] = InZriH, J] (120) 

where 

N N 

i=i -'^ • 

and 

Zjv[//, J] = rr^e^-^-*°--^-^ . (122) 
Functional derivatives with respect to J generates a factor of g in an average 

{ga) = ^W[H,J] (123) 

C. Perturbation Theory and Hierarchial Structure of Theory 

It is trivial to see that the fundamental identity of the theory takes the form 
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G„ = TrVae^-^o+-^%^^[^''^l (124) 

where 

^W[H-J] = W[H + F-J]-W[H-J] (125) 

in the extended space. 

We can again construct a dynamical OZ equation of the form 

Gap = QaP + Cafj.G^i3 (126) 

where the single-particle contribution is given by 

Qap = TrV,0^e^-<^°+^-%^^[^^^l (127) 
and the dynamic direct correlation function 

= TrVae^-*°+^ ^'e^^[^'^]-^AW^[/7; J] (128) 
The key point is that to first order in perturbation theory 

AW[H;jf^ = J2F^,Guo (129) 

Uo 

where the sum over uq is only over the core variables 0o = {p,B}. This happens because 
the Hamiltonian and hence the corresponding MSR action involves only these core variables. 
The first order dynamical direct correlation function c^^^ vanishes ii fi = g. This means at 
first order we have the kinetic equation 

Gal3 = Qal3 + Gapl^VGpp + GaBPVGBp ■ (130) 

The structure of this equation is very interesting. If we restrict a and (3 to the core variables 
then Eq. fll30p reduces to the first- order core problem we solved earlier. The correlation 
functions for the extended variables $ can be expressed in terms of (roughly) noninteracting 
correlation functions between all the variables and solutions to the core problem associated 
with $0- 



25 



This suggests that one can not get around the core problem. This problem is self- 
contained and must be treated first. These self-consistent solutions then enter into the 
extended space $ as determined sub-matrices in the matrix kinetic equation. The second 
order order corrections to the dynamic direct correlation functions c\j and hence the self 

(2) 

energy matrix K^, will involve the extended set of collective variables. 



VI. NON-INTERACTING SYSTEM 



A. Generating functional 

In this section we present the calculation of the zeroth- order cumulants for the Newtonian 
dynamics case. We first work out the generating functional for a rather general single- 
particle gaussian model in the Appendix [Bl This model includes SD, ND and Fokker-Planck 
dynamics as special cases. For a quadratic action of fields ipi and response fields ipi we have 



^ rOO ^ _ ^ /"OO 

r -s ^ 

- J2 h,{t)^i{t) + hi{t)ij,{t) 



(131) 

where Dij is the damping matrix, Ky is a force matrix, and and hi are the detailed 
external fields that couple to ipi{t) and ipi(t). In Appendix [Bl we find that the associated 
generating functional is given by 



lnZo(/i,/i;^(°)) = J2jd^J dt'\hi{t)c^jit-t')hj{t') + hi{t)gij{t-t%{t') 
ij 

+ dthi{t)igij{t-to)i;^ 



',(0) 



(132) 



where ipf'^ is the initial value of the fields, the normalization is such that Zq[0, 0] = 1, 



/oo _ _ _ 

dtgik{t,i)Dk,gJjit,t') 
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/oo _ _ _ 

dig,e{t', t)2DugUi, t') = c,,{t', t) (133) 

and the function gij is now obtained from the solution of the Green's function equation 

d 

g^dijit, t') + ^ KikQkjit, t') = -i6{t - t')6ij . (134) 

These equations govern SD, ND and FPD. Notice that all information about the equilibrium 
state of the system is carried by the damping matrix D ^ ksT. If we restrict ourselves to 
simple fluids, then the damping term with Cij is zero and the state of the system enters via 
the initial condition. We have 



P OO /"OO 

Zo{h,h,tlj^^^) = exp{ dt dt'J2Ht)9ij{t-t')[hj) + i6{t' -to)ijf]) (135) 
with an average over initial conditions remaining. If the initial conditions are gaussian 



P^^ifi))=Ue-^^^^^^' (136) 



one can carry out the average over ijj^^^ and obtain 



^ 1 ^ 
InZoih, h) = - dt dt'Y,hi{t)cl{t - t')hj{t') + hi{t)gij{t - t')hj (137) 

2 J to J to - 



and 



Clit,t') = Q,(t,t')+4(^'^') 

4{t,t') = Y.9^kit,t,){M-%,gl{to,t') . (138) 
B. Newtonian Dynamics 

For Newtonian dynamics there are two phase-space coordinates, R and P, with an ele- 
mentary force matrix Kij = 5i^R5j^p^ for the noninteracting system. We first have to solve 
the matrix equation following from Eq.f ll34p in the Newtonian case for which indices i and 
j run over the set {R, P}. 
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R 


P 


R 


-iOit - 1') 




P 





-i9{t-t') 



TABLE I: Matrix g: Newtonian Case. 



d 


9pr 


Q^gRR 


m 


d 


9pp 


^^^^ 


m 




d 




^^^^ 




d 







= -i6{t) (139) 

= 

= 

= -iS{t) (140) 



The straight forward solution gij{t — t') is hsted in Table 1. 



Since there is no dissipation put into the model, D^j — 0, cf^ — the generating functional 
simplifies 



ZnZo[/i,/i;V'f] = / dt dt'Y,hi{t)gij{t-t')hj{t') 

Jto Jto 

oo 

hR{t)e{t - to)Ro + % E hi{t)i9ip{t - to)Po 

i 

The general result for the solution of the differential equation 



+ [ dt 

Jto 



(141) 



d_ 
di 



(142) 



is obtained in terms of the Green's function g{t, t') 
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< ^i{t; to) >= E ^ / dt'9^,{t - t')i^T5{t' - to) = ^g^J{t - to)^f • (143) 



For Cij = and h = h = we obtain the same relation from eqn. f lB12p . In this particular 
case, we obtain the following reversible ( and deterministic) equations of motion are given 
for the two coordinates {R,P}. 



R(t) = igRR(t-to)Ro + igRp(t-to)Po 
P{t) = igpR(t-to)Ro + igpp(t-to)Po 

These equations are also obtained from the relations 



(144) 



< R{t) > 

< m > 



6hR 

6 
6hp 



InZo 
InZo 



h,h=0 



h,h=0 



Inserting the g^s leads to the standard phase-space trajectories for free streaming particles. 



R(t) = 9{t-to 
P{t) = e{t~to)Po 

The generating function for ND is given by 



Ro + it-to) — 
m. 



(145) 



lnZo[h, h] 



oo pca 



dt / dt'J2hi{t)gij{tt')hj{t') 



roo roo 

+ / dt{hR{t)e{t-to)Ro+ / dtY,K{t)gip{t-to)Po . (146) 

J to J to ■ 



C. Computation of ^-Correlations 

It is not the direct correlations of the phase-space coordinates that are of interest but 
the correlations of the $. It is clear that we need to evaluate all of the p — B cumulants 
generated by 



Wa[H] =Tre 



(147) 
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We begin by introducing the microscopic sources hi, hi and treating 

Zo[H,h,h]^Tre"-'l'e''-'f'+'^-'f' . (148) 

Remember that the spatial Fourier transforms of the are in terms of the phase-space 
coordinates ■0: 



0p(l) = e-*-^(*i) 

(Pb{1) ^ -ki ■ P{ti)e-''''-^^''^ . (149) 

Introducing the operators 

0^(1) = e"*^'^ (150) 
0b(1) = 6(l)<^p(l) (151) 

(152) 

where 

6(1) = -^1--^ • (153) 

ohp{ti) 

We can write 

Zo[H, h, h] = rre^-^+'^-'^+'*-^ (154) 
= e"-hre^-'l'+'^-^ (155) 
^e^-'^poe^-^-^^ j djPoP[iPo]e''-'''''° (156) 

= e^'Vo / ciV^oPl^ole'-^-t^+^^ol (157) 

One can then express an arbitrary noninter acting cumulant, in the presence of h and h, 
in the form 

Gbb...Bp..M^ ...U+l...n;h,h)= 6(1) . . . 6(£)0(1) . . . 0(n)Zo[/i, h] (158) 
= 6(1) . . . 6(£)0(1) . . . 0(n)po / d2PoP[He''-'-^^^''^°^ (159) 
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The analysis is very similar to that for the SD case in Appendix C in Ref. I4]. The 0p are 
translation operators: 



GBB...Bp...p{l2 ...a + l...n;h,h) = b{l) . . . b{i)po I (i^oP[V'o]e["+^l-^-[^+''^°l (160) 
where in more detail 

h + L ^ K{t,) + L^^\t^) (161) 

with 

L^a\t,) = -i5^,R jz ks5{t, - t.) (162) 

s=l 

Application of the operators h{i) is multiplicative 

GBB...Bp...p{l2..Ai+l...n-h,h)=hn{l)..M(^)po J rfV^oPWef'+^l-^t^+^^^l (163) 

and 

roo _ _ _ 

bnU) = -kj / diJ2 [ha{i + Ll{t)] g^p{i, tj) . (164) 

J to Q, 

Thus we have the cumulants for arbitrary h and h. Setting h and h to zero we have 

GBB...Bp...p{^2...U+l...n)=b^{l)...bn{i) .po J d^oP[He''''"^°^ (165) 
= 6„(1) . . . bn{i)po J d^PoPiHe^"^' (166) 

where 

K^ = iL-g = tL^;\i)g^^{t,to) (167) 

and 

bn{j) = -k, / dii-iJ2kJii-Q)i 9it-tj)ii-t,)) (168) 

k ■ " 

= - J2''s9its-tj){ts-tj) (169) 
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We have after a little work 



n 



«i^ = -^E^3 (170) 



s=l 



Kp^-iJ2^s{ts-to) (171) 



and 



= -iT (172) 

Gbb...Bp..A^2 ...ii + l...n) = bn{l) . . . hn{t) (173) 

X po/ dfiRQP[RQ\e'^"=^^^-^° J d'^Po^[^o]e~'^-^° (174) 

= bn{l) . . . hn{€)po{'2nY k)j e^" (175) 



where 



^ J d'^PoPiPole-'^-^^ . (176) 

The probability of finding a particle at a given position is uniform and the position average 
enforces the translational invariance in the system 

I ci'='i?oe^«-^° = {2n)'^5{h + k2 + ... + K) . (177) 



The momentum is assumed to obey a Maxwell-Boltzmann distribution and the average gives 
the contribution 

J (27rpo)'^/' ^ ^ 

where Pq = 'mf3~^ — iti^Vq is the thermal momentum. This is a complete the square calcu- 
lation giving 



e^" = e-^^o^ (179) 

and we identify the effective 

Nn = -lplT' . (180) 
The correlation function Gbb...Bp...p with i number of B fields is obtained as 
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GBB...Bp...p{l2 ...££ + l...n) = bn{l) . . . bn{i)po{2rrY [J2 ^sj e""" • (181) 
Let us focus on A^„. Writing it out we have 



2m? 



^ ks{ts - to) 



ls=l 



(182) 



Using the conservation law Eq. (I177p we see this expression is independent of to and is time 
translationally invariant. The effective can be written in a form similar to the SD case. 
We can rewrite Nn 



where 



Then define 



We then have 



Nn — — - X! X! ^i^j^ij 
^ i=l j=l 



r - ^" f f 



v',it,-t,y 



(183) 



(184) 



(185) 
(186) 



^ /c^Cji + ^ ^ kikjCij 

i=l i=^j=lj=l 
n "1 

+ X! kikj-[Cii + Cjj — Dij] 
i=i i^j=i ^ 

n n n "1 

^^k'^Cii+ ^ kikj-Cii+ ^ kikj-[Cjj+ ^ kikj-{—Dijj 

1=1 i¥^j = ^ «7^i = l *7^i = i 







2 

2 n rt 



n "1 "1 

^ fcj Cji — 2 ^ /c- -Cjj + ^ kikj-{—Dij] 

i=l i=l ^ i+3=^ 



^ i=\j=\ 



(187) 
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where 

Kij = ki ■ kj . (188) 

The noninteracting correlation and response cumulants are now all in a time-translationally 
invariant form. 

D. Two- Point Cumulants 

Let us extract the two-point cumulants needed in our first-order discussion. For the 
density-density correlation function 

Gf}{l2) = p^{2T:Y5{k^ + A;2)e-^'^^"o(*i-*^)' (189) 

where 

Ki^kl . (190) 

For the response function we have 

G%{l2) = h,{2)Gf^{l2) (191) 

where 

62(2) = --^(ti-t2)(ti-t2) (192) 
m 

and 

4^(12) = 62(1)62(2)G(°)(12) = . (193) 

Notice the relationship 

G%{l2) = e{U-t,)P^Gf^{l2) . (194) 
holds at the zeroth order. We return to these results in the next paper in this series. 

VII. CONCLUSIONS 

We have developed the fundamental theory for conventional fluids in which the particles 
follow reversible Newtonian dynamics. The theory is remarkably similar to that developed 
to treat particles following Smoluchowski dynamics(SD). We have a field theoretic formu- 
lation with similar nonlinear interactions in terms of a pair potential. The question of 
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self-consistency is addressed in a very similar manner. The role of initial conditions and 
broken time translational symmetry are addressed at zeroth and first order in the theory. 
At zeroth order the average over initial conditions is shown for cumulants of the collective 
variables to satisfy TTI. At first order wc can impose equilibrium and sustain TTI if we 
require the system to obey fluctuation dissipation symmetry (FDS). 

As for the treatment of SD the second order self-energies require the zeroth order irre- 
ducible three-point vertices. These three-point functions turn out to be more complicated 
in the case of ND. In the present paper we have shown how to compute all cumulants for 
the non interacting system. The three-point cumulants enter into the determination of the 
three-point irreducible vertices. These vertices enter into the computation of the second or- 
der self-energy. This becomes the necessary input for analysis of a possible ENE transition 
in such systems at high density. In the next paper in this series we show that there is a 
FDS and it can rather easily be applied to the cumulants and irreducible vertices of the 
fully interactive fluid ND system. We derive identities obeyed by the three-point objects 
in a manner similar to the application of nonperturbative FDS to the two-point cumulants 
and vertices. We then focus on the use of this machinery to explore whether we have ENE 
transitions in ND systems. 
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Appendix A: Fundamental Identities 



We have the definition of the grand partition function 



oo 



Zt[H;x] = Yl :^rr^e^-*^+^*^-*^ (Al) 

where — J2iLi4'^^^ and we have the important result that self-energies ( the interaction 
of the i-th particle with itself does not contribute, i.e., 

= (A2) 

for all i. We need the following identity 

$JV+I(7$jv+1 = ^N(T^N + 2(/)(^+^)(7$jv = ^N(T^N + 2F^+^$jv (A3) 

where we have defined F^^'^^'^ = acp^^^^K Therefore we have 

^^N+ia^N+i = ^*jv(j*jv + F(^+i)$jv ■ (A4) 

Now consider 

Zj,{H + F^+i) = Tr^e(«+^'"+'')-*-+^*-'^*- (A5) 
Multiplying by e^''^'^+^ gives the result 

This is the canonical form of the central identity. We now Trace over the degrees of freedom 
of the {N + l)-th particle ( which we denote by the label 0, i.e., N ^ 0) 



= Zm+i{H) (A8) 
where Tr^'^^ implies the trace taken over the phase space coordinates of the 'O'-th particle. 

N 

We now multiply the above equation by ^ and sum over all N 
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oo N oo ™JV-1 

N=0 ■ N=l 

Divide the above equation by Zt[H] x] to obtain 

^W[H; x] = fr(o)^H.^(o,+AW[H+F(oy.] (^^^q) 
where In Zt[H] x] = W[H; x] and we have defined 

AW[H + F;x] = W[H + F^°^;x]-W[H;x] . (All) 
Integrating the relation ( ]A9I) we obtain 

l^(/7;po)= rrfxTr(°)e^-^™+^^[^+^^°'^^] . (A12) 
Jo 

We use this result to generate the one-point quantity 
6W S fpo 



Jo 

+ r rfa;f rWe^-^'°'+^^(^+^'"'^^) iG^iH + F^; x) - G^{H- x)] . (A13) 

Jo L -I 

Substituting H' = H + F^^^^ the integral in the second term on the RHS is written as 

r rf2;f r(0)e^-^*"'+^^(^ [G^iH'; x) - G^{H; x)] . (AM) 
Jo 

Consider first the first part as, 



h = r dxfT^'^e'' 'f'"'\'''^^''''^^G^{H'- 
Jo 

Zt{H';x) 
Zt{H;x) 



x] 



r dxfr^%^'^'^^G^{H'- X) (A15) 



We want to show self-consistently that 



37 



satisfies the eqn. (]A13P above. Substituting eqn. (1A16P in eqn. (lAlSP we obtain 



Zt[H) 



where in getting the last equahty we have used the result 



Using the last relation (lA9p . we obtain from eqn. (1A17P 



0« 



The expression for I2 is put back into the expression (lAlSp for Ga- 



Jo 

po 



Zt[H + F^^^] 
Zt[H;x] 



Jo 

= Tr(^)0«e^-*'^' r dx 



JO ox [ Zt[H; x\ 
Zt[H;po] 



ZT[H-x]dx ^ Z^[H;x] dx ^ 

Zt[H + F^'^'>] d Zt[H + F^^^] 

Zt[H] Zt[H;x] 
Zt[H + F^^^-x\ 
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Thus we have estabhshed the fundamental result 

^ j^^il)^il)^H-<t>W+AWlH+Fm] 



(A20) 



Appendix B: Single-particle Gaussian Problem 

We want to construct the generating functional associated with single-particle noninter- 
acting dynamics. The single-particle problem is governed by the phase-space coordinates ipi 
with response variables The action governing these variables is quadratic 



to 



dt 



to 



?1 



dt 



hi{t)jPi{t) + hi{t)i^i{t) 



(Bl) 



where Dij is the damping matrix, Kij is a force matrix, and hi{t) and hi{t) are the detailed 
external fields that couple to particle ip and ip. Newtonian dynamics corresponds to the 
special case where the damping matrix vanishes Dij — 0. 

We proceed using the identities that hold in the range to < t < oo: 



J 6ipi[t) 

which leads to the set of equations where we sum over repeated indices labeled by j and we 
suppress the local time label: 

2DijGj + iGi + iKijGj = h, 
-i^Gi + iKfjGj = hi 

where is the transpose of the matrix K. We have defined the functions Gj and Gi as 



Gi 
Gi 
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The Gs depend on initial data which we must eventually average over. We must now solve 
these equations, at least formally, to obtain the generating functional. 

Let us treat Gi(t) and Gi(t) as inner products spanned by a complete and orthonormal 
set of states \i > such that 



G,{t) 


= < i\G{t) > 


G,{t) 


= < t\G{t) > 




= < i\K\j > 




= < i\D\j > 


hi{t) 


= < i\h{t) > 


kit) 


= < i\h{t) > 



which introduces the operators and vectors G,G, K, D, h, h. We then have the operator 
equations 



2DG{t)+t(^^^+K^ Git) 



hit) 
hit) 



We first solve the eqn. IB3I for G, obtaining 



Git) 



drg' (t,r)/i(r) 



(B2) 
(B3) 



(B4) 



where 



g^it.t') = -ie-'^''eit' -t)e 



+ K 



Git) = -i (hit) - 2D p dt'g^it,t')hit'] 



The solution for G is put back in the eqn. flB2l) for G to obtain 

d_ 

dt 

It is straightforward to obtain the corresponding solution for Git) as 

Git) = igit, to)V^^°^ + / dt'igit, t')hit') + c(t, t')hit')] 



(B5) 



(B6) 



(B7) 
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where the functions g and c are given by 



g{t,t') = -ie-^*0(t - t')e^*' (B8) 

/oo _ _ _ 

dtg{t,t}2Dg^{t,t') (B9) 

and ^/'*^°^ is the initial value of the phase-space coordinates. Notice that g satisfies the Green's 
function equation 



g{t,t') = -i5{t-t') 
Putting in complete sets of states we obtain 



(BIO) 



G^{t) = Y.j dt'h,{t')g,Xt' -t) (Bll) 
i 

G^(i) = T.r ^^'^3^,it - t')kit') + Cy(t - t')h^{t')] + ig,,{t - to)^ljf (B12) 



where the kernel matrix Cjn is obtained as 



/oo _ _ _ 

dtgik{t,t)2Dk,gJj{t,t') 

/OO _ _ _ 

dtgje{t',i)2DM9li{t,t') = c^iit' ,t) 
. ., -oo 



(B13) 



In getting the last equality we have used the symmetry Dij = Dji of the damping matrix. 



The function gij is now obtained from the solution of of eqn. (IBlOp with labels restored 



d_ 

di 



9ij{t, t') + K,kgkj{t, t') = -i5{t - t')5i_ 



We then have the results for the generating functional 

5 



In Zo(^) h) 



5h(t) 



(B14) 

(B15) 
(B16) 



The solution to this set of equations, Eqs. (]Bll|) . (]B12p . (!29|) and (lB16P for the generating 
functional is given by 
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+ dthi{t)igij{t-to)^l:f 

= ]-h- c-h + h- g -h + h-ig ■ . 



(BIT) 



The last short-hand representation is useful. The full generator requires averaging over the 
initial conditions 



exp 



1 



- h • c • h + h ■ g ■ h + h ■ ig ■ if] 



(0) 



= 62 



hh-c-h-\-h-g 



(B18) 



This is the solution to a rather general gaussian problem. 
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